Determining Equation of State Binary Interaction Parameters Using K- and L- Points

In presenting this thesis in partial fulfillment of the requirements for a Postgraduate degree from the University of Saskatchewan, I agree that the Libraries of this University may make it freely available for inspection. I further agree that permission for copying of this thesis in any manner, in whole or in part, for scholarly purposes may be granted by the professor or professors who supervised my thesis work or, in their absence, by the Head of the Department or the Dean of the College in which my thesis work was done. It is understood that any copying, publication, or use of this thesis or parts thereof for financial gain shall not be allowed without my written permission. It is also understood that due recognition shall be given to me and to the University of Saskatchewan in any scholarly use which may be made of any material in my thesis. ABSTRACT The knowledge of the phase behaviour of heavy oils and bitumen is important in order to understand the phenomenon of coke formation. Computation of their phase behaviour, using an equation of state, faces problems due to their complex composition. Hence n-alkane binaries of polyaromatic hydrocarbons are used to approximate the phase behaviour of heavy oils and bitumen. Appropriate values of binary interaction parameters are required for an equation of state to predict the correct phase behaviour of these model binary fluids. This thesis deals with fitting of the binary interaction parameter for the Peng-Robinson equation of state using landmarks in the binary phase space such as K-and L-points. A K-or an L-point is a point in the phase space where two phases become critical in the presence of another phase in equilibrium. An algorithm to calculate K-and L-points using an equation of state was developed. The variation of calculated K-and L-points with respect to the binary interaction parameter was studied and the results were compared with the experimental data in the literature. The interaction parameter was then fitted using the best match of experimental results with the computed ones. The binary interaction parameter fitted using a K-or an L-point was then used to predict the P-T projection of the binary system in phase space. Also, the qualitative effect of the binary interaction parameter on the P-T projection was studied. A numerical and thermodynamic study of the algorithm was done. Numerical issues like the initial guesses, convergence criterion and numerical techniques …